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ABSTBACT 

Eecent  improvements  in  rare  earth  magnets  have  made  it 
possible  to  construct  strong,  lightweight,  high  horsepower 
EC  motors.  This  has  occasioned  a  reassessment  of  electrome- 
chanical actuators  as  alternatives  to  comparable  pneumatic 
and  hydraulic  systems  for  use  as  flight  control  actuators 
for  tactical  missiles.  A  dynamic  equivalent  circuit  model 
for  the  analysis  of  a  small  four  pole  brushless  DC  motcr  fed 
by  a  transistorized  power  conditioner  utilizing  high  speed 
switching  power  transistors  as  final  elements  is  presented. 
The  influence  of  electronic  commutation  on  instantaneous 
dynamic  motor  performance  is  particularly  demonstrated  and 
good  correlation  between  computer  simulation  and  typical 
experimentally  obtained  performance  data  is  achieved.  The 
model  is  implemented  in  CSMP  language  and  features 
improvements  in  transistor  and  diode  models  as  well  as  a 
more  accurate  air  gap  flux  representation  over  previous 
work.  Eall  effect  sensor  rotor  position  feedback  is  simu- 
lated. Both  constant  and  variable  air  gap  flux  is  modeled, 
and  the  variable  flux  model  treats  the  flux  as  a  fundamental 
and  one  harmonic. 
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I.  INTBODOCTION 

Electronically  ccmmutated  brushless  DC  motors  are 
increasingly  becoming  practical  for  numerous  applications. 
The  advantages  of  these  motors  over  conventional  brush-type 
IC  motors  are  smaller  size  and  lighter  weight  per  horse- 
power. The  increase  in  power  over  conventional  DC  motors  is 
realized  by  the  use  of  three  phase  circuitry  commonly  found 
in  synchronous  AC  motors  which  results  in  additional  motor 
current-magnetic  flux  interaction  not  available  with  conven- 
tional DC  motor  designs.  The  elimination  of  the  trush 
commutation  and  rotating  armature  allows  operation  at  higher 
speed  with  higher  currents  and  improves  the  thermal  charac- 
teristics of  the  motor  as  heat  generating  windings  will  be 
contained  in  the  stator  which  can  be  more  efficiently 
cooled. 

In  this  thesis,  an  equivalent  circuit  motor-power  condi- 
tioner modeling  approach  is  used  to  predict  the  dynamic 
performance  of  a  typical  brushless  DC  motor  for  a  fin 
control  actuator  assembly.  The  interaction  of  the  power 
conditioner  final  drive  elements  with  the  stator  currents  is 
demonstrated  in  addition  to  the  rotor  magnetic  flux  -  stator 
current  interaction  and  interdependence.  Numerical  results 
of  steady  state  operation  are  compared  with  experimentally 
obtained  instantaneous  voltages  and  currents  in  amplitude 
and  in  profile.  After  demonstration  of  steady  state  opera- 
tion of  the  model,  start-up  and  reversal  transients  are 
predicted. 

The  model  developed  in  this  thesis  is  part  of  a  contin- 
uing program  to  develop  a  comprehensive  simulation  for  a  fin 
control  actuating  system  which  has  been  the  subject  of 
previous  work  [Eef.  1  and  2]  and  is  subject  to  further 
development. 
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The  object  of  this  thesis  is  to  extend  previous  work  on 
brushless  DC  motors  [Bef-  1 ]  to  include  switching  character- 
istics of  power  transistors  and  diodes  in  such  a  manner  as 
to  provide  a  coherent  model  of  power  conditioner  final  drive 
elements  and  a  brushless  DC  motor.  Additionally,  the  simu- 
lation provides  steady  state  as  well  as  instantaneous 
dynamic  performance  during  start-up  and  sudden  reversal  of 
the  irotor.  The  IEH  Continuous  System  Modeling  Program 
(CSMP)  was  chosen  primarily  because  of  the  inherent  attri- 
butes of  CSMP  for  modeling  dynamic  physical  systems. 
Inasmuch  as  previous  model  development  utilized  this 
language,  seme  continuity  of  development  is  maintained. 


II.  SYSTEM  DESCRIPTION 

A.   SISTEH  BLOCK  DIAGBAH 

The  simulation  diagram  for  the  model  is  depicted  in 
Figure  2.1  and  demonstrates  the  interrelation  between  the 
magnetic  flux  model,  the  torque  coefficient,  and  the  gener- 
ated tack  eiif  voltage.  The  nonlinear  nature  of  the  gener- 
ated voltage  and  torque  is  indicated  by  the  multipliers  in 
the  diagram  for  angular  speed  times  the  magnetic  flux, 
M  (t)  (J)  (t)  ,  used  in  the  back  emf  calculation  and  the  motor 
current  times  the  magnetic  flux,  I(t)Q(t),  used  in  the 
generated  torque  calculation.  The  three  phase  electrical 
simulation  calculates  three  phase  currents,  which  are  multi- 
plied by  the  torque  coefficient  and  normalized  flux  nodel 
for  each  phase  resulting  in  a  torque  due  to  the  current  in 
each  phase.  The  motcr  torque  is  the  sum  of  the  three  phase 
torques,  and  the  net  torque  is  the  motor  torque  less  the 
load  torque.  Motor  rotation  is  then  calculated  by  applying 
the  net  torque  to  the  rotor  inertia  and  viscous  friction. 
The  integral  of  the  rotor  angular  speed  is  calculated  for 
angular  position,  which  is  fed  to  the  flux  and  transistor 
logic  models  to  control  the  generation  of  conduction 
currents  for  each  phase.  The  transistor  model  and  magnetic 
flux  model  are  discussed  in  Chapter  Three  and  Chapter  Four, 
respectively. 
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TABLE  I 

Typical 

Motor  Parameters 

Parameter 

Value 

Units 

Stator  resistance,  Ea 

1.37 

Ohms 

Stator  Inductance,  La 

0.0016 

Henrys 

Torgue  coefficient,  Kt     15.9 

Oz-In/amp 

Eeverse  voltage 

coefficient,  Kb 

0.112 

Volt-Sec/Ead 

Ector  inertia,  Jm 

0.001 

Oz-In/Sec 

Viscous  friction 

coefficient,  Em 

0.022 

Oz-In-Sec/Ead 

L 

.,  ._ 

B.   MCTOE  CONFIGURATICH 

The  motor  simulated  in  this  thesis  is  a  three  phase, 
four  pcle  device  typical  of  modern  brushless  DC  motors  as 
described  in  [ Bef .  1  ].  Typical  values  for  the  motor  parame- 
ters are  contained  in  Table  I.  The  listed  values  of  resis- 
tance and  inductance  are  as  measured  across  two  windings  as 
the  motor  is  only  a  three  wire  device  as  pictured  in  Figure 
2.2.  It  should  be  noted  that  the  motor  inertia  is  very  low; 
a  typical  rotor  diameter  is  of  the  order  of  one  inch  and 
results  in  a  very  small  mechanical  time  constant.  A  tran- 
sverse cross  sectional  view  of  rotor  and  stator  for  a  the 
model  motor  is  shown  in  Figure  2.3,  with  the  rotor  in  the 
zero  (or  180)  degree  position.  The  magnetic  rotor  can  be 
seen  to  have  north  and  south  poles  separated  by  ninety 
degrees.  This  generates  a  flux  distribution  which  alternates 
polarity  every  ninety   mechanical   degrees   as  depicted  in 
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Figure  2,2    Hctor  and  Ccmmntating  Transistors 

Figure  4.4.  Since  tie  phase  windings  are  separated  ty  sixty 
mechanical  degrees,  the  flux  in  each  phase  is  seen  to  have  a 
sixty  degree  phase  shift  with  respect  to  the  adjacent 
phases.  The  explanation  of  the  shape  of  the  model  flux  is 
discussed  in  Chapter  lour.  The  electronic  commutation  must 
he  performed  in  such  a  manner  as  to  complete  two  electrical 
cycles  for  each  mechanical  rotation  since  the  flux  distritu- 
tion  completes  two  cycles  in  one  rotor  rotation.  The  commu- 
tation of  the  armature  (stator)  current  is  accomplished 
every  thirty  degrees  of  rotor  rotation  because  a  new  winding 
"enters"  the  magnetic  field  with  every  thirty  degrees  of 
rotation  simultaneously  with  another  winding  "leaving"  the 
magnetic  field-    The   sensor  and  switching  logic   for  motor 
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Figure  2.3    Motor  Stator  and  Winding  Configuration 

rotation  is  tabulated  in  Table  II,  and  indicates  the  sensor 
logic  level  (EPS),  and  connection  of  the  phase  ternirals  for 
both  clockwise  and  counterclockwise  direction  for  each 
thirty  degrees  of  rotor  position.  Note  that  regardless  of 
direction  of  rotation,  the  rotor  position  logic  remains 
constant  for  a  given  rotor  position,  and  that  rotor  posi- 
tions between  180  degrees  and  360  degrees  are  indistingui- 
shable from  rotor  positions  tetween  zero  and  180  degrees 
because  of  the  rotor  symmetry. 
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TABLE  II 

Sensor  and  Switching  Logic 

ECTOE 
POSITION 

EPS 
A 

Cou 

BES 
B 

nterclockwise  flotation 

EPS     PHASE      PHASE 
C        A          B 

PHASE 
C 

0 

30 

60 

90 

120 

150 

0 
1 

1 
1 
0 
0 

1 
1 

0 

0 
0 

1 

0      open      neg 
0      pos       neg 

0  pos       open 
open      pos 
neg       pos 

1  neg       open 

pes 

open 

neg 

neg 

open 

pos 

EOTCE 
POSITION 

EPS 
A 

EPS 
| 

Clockwise  flotation 
EPS     PHASE      PHASE 
C        A           B 

PHASE 
C 

0 

3C 

60 

90 

120 

150 

0 
1 
1 
1 
0 
0 

1 
1 

0 

0 
0 

1 

0      open      pos 
0      neg       pos 

0  neg       open 
open      neg 
pos       neg 

1  pos       open 

neg 

open 

pos 

pos 

open 

neg 

C.   HCTOB  SYSTEM  EQUATIONS 

Derivation  of  tie  basic  motor  system  equations  (equa- 
tions 2.1  through  2.8)  was  accomplished  by  Thomas  in 
[Eef.  1]  and  they  are  repeated  here  for  convenience.  From 
Kirchoff's  voltage  law,  the  stator  voltage  must  be: 


Es(t)  =  Ldl/dt  +  EI(t)  +  Eb(t) 


(2.1) 


where  the  back  emf  voltage  can  be  described  by: 


Eb  (t)  =  (Km  |)flm(t)  =  KbWm(t) 
or   Eb(t)  =  Kbd«m(t)/dt 
since   dGm(t)/dt  =  Wm  (t)  . 


(2.2) 
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Rearranging  Equation  2. 1  for  the  highest  derivative  on  the 
left    side   yields: 

dl/dt    =     (1/L)Es(t)-(B/L)I  (t) -(1/L)  KbWm  (t)  .  (2.3) 

The   mechanical   torque  balance   equation  is: 

Tm  (t)    =   JdWm(t)/dt   +   Bflm(t)    ♦  Tl(t).  (2-4) 

Rearranging  yields: 

dWm(t)/dt  =    (1/J)Tm(t)    -    (1/J)Tl(t)    -    (B/J)  Wm  (t)  (2.5) 

The  total  viscous  friction  (B)  and  total  load  inertia  (J)  is 
represented   as   follows: 

B    =    Bin    ♦    El  (2.6) 

J  =  Jm  +  Jl  (2.7) 

where  Bl  =  Blp/U  and  Jl  =  Jlp/N.  Blp  is  the  viscous 
friction  associated  with  the  load,  Jlp  is  the  load  inertia 
and  N  is  the  gear  reduction  ratio.  The  link  between  the 
electrical  and  mechanical  equations  is  the  torque  constant, 
which  is  derived  from  the  magnetic  flux: 

Tm  (t)  =  Km$I  (t)  =  Ktl  (t)  .  (2.8) 

The  basic  equations  have  been  enhanced  in  this  model  to 
take  advantage  of  inherent  sign  conventions,  particularly  in 
direction  of  current  flow.  This  in  turn  makes  the  commuta- 
tion procedure  more  apparent  in  the  generated  graphics.  lor 
example,  the  primary  motor  currents  are  depicted  in  Figure 
2.2,  each  having  an  assigned  positive  direction  as 
pictured.   The  difference  between  TAB  and  ICA  yields  IA,  IBC 
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minus  IAE  equals  IB  and  ICA  minus  IBC  equals  IC.  The  current 
flow  cannot  always  be  positive  as  depicted,  since  the  sum  of 
the  three  phase  currents  must  be  identically  zero  by 
Kirchcff's  current  law.  Therefore,  at  least  one  of  the 
currents  IA,  IB  or  IC  must  be  negative  at  any  given  time, 
representing  a  flew  of  current  returning  to  the  power 
supply.  From  the  current  waveform  it  is  immediately 
apparent  which  phase  is  connected  to  the  negative  supply  and 
which  is  connected  to  the  positive  supply  simply  by 
observing  the  polarity  of  the  current  (see  Figure  5.2).  The 
concept  of  the  positive  and  negative  supply  voltages  is 
totally  eguivalent  tc  a  single  positive  (or  negative)  supply 
as  used  in  [Eef.  1]  and  since  the  drive  currents  are  calcu- 
lated from  the  difference  of  two  fictitious  currents,  a 
single  supply  will  work  in  the  simulation,  but  the  split 
supply  enhances  the  symmetry  of  the  motor  currents  and  volt- 
ages. This  change  allows  for  improved  modeling  and  conceptu- 
alization of  transient  studies  as  discussed  in  the  next 
secticn. 

E.   HCIII  DEVEIOPMEHT 

The  tasic  simulation  of  the  brushless  DC  motor  by  Thomas 
used  a  single  positive  or  negative  power  supply  and  superim- 
posed armature  leg  currents  to  produce  the  measured  motor 
performance.  Motor  drive  was  then  realized  by  multiplying 
averaged  armature  currents  by  a  torgue  factor  for  the  model 
under  development.  In  the  this  model  the  supply  is  a  split 
supply  of  equal  voltages.  Armature  current  is  assigned  a 
positive  sign  if  it  flows  in  the  positive  direction  (i.e. 
into  the  motor)  as  shown  in  Figure  2.2  and  a  negative  sign 
if  it  flows  in  the  opposite  direction  (i.e.  out  of  the 
motor)  as  discussed  previously.  A  secondary  result  of  the 
model  of   commutating  transistors  and   split  supply   is  that 


17 


motor  reversal  cannot  be  accomplished  by  reversing  the 
polarity  of  the  power  supply  as  in  £Be£.  1]  but  must  instead 
be  accomplished  as  in  a  typical  power  conditioner  by  proper 
reverse  sequencing  of  the  commutation,  shown  in  Table  II, 
and  discussed  in  Chatter  5. 

Several  models  were  produced  whereby  each  model  repre- 
sented a  completion  cf  a  particular  aspect  of  the  system 
such  as  motor  air  gap  flux  and  power  conditioner  transistor 
switching  characteristics.  The  early  models  were  developed 
from  simple  ideas  such  as  instantaneously  switching,  zero  or 
infinite  resistance  transistor  switches  and  constant  air  gap 
flux.  Each  model  improved  a  particular  area  under  develop- 
ment until  it  was  thought  to  be  a  suitable  representation  of 
an  actual  motor  or  power  conditioner  feature. 

The  first  model  incorporates  a  square  wave  of  back  emf 
(constant  flux)  and  endeavors  to  demonstrate  transistor 
commutation  effects.  The  transistor  switching  is  instanta- 
neous and  is  accomplished  within  the  "commutation  and  Hall 
effect"  procedure   (see  Appendix  C)  . 

The  second  model  introduces  changes  in  the  transistor 
switching.  While  switching  is  still  instantaneous,  typical 
values  of  transistor  saturation  and  cutoff  resistance  are 
introduced,  and  commutation  is  controlled  by  the  commutation 
procedure  instead  of  being  accomplished  by  that  procedure. 
This  feature  leads  to  the  next  model  which  simulates  the 
actual  rise  time  and  saturation  delay  of  power  transistor 
switching.  This  model  is  the  basis  for  all  subsequent  simu- 
lations involving  increased  model  complexities  such  as  the 
sinusoidal  flux  model. 

The  remaining  simulations  develop  the  voltage  and 
current  relationships  between  the  motor  and  electronic 
commutator  final  drive  transistors  and  incorporate  sinu- 
soidal magnetic  flux  nodels  consisting  of  a  fundamental  and 
fifth  harmonic  as  described  in   Chapter  Four.   These  program 
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listings  are  contained  in  Appendix  B.  The  first  of  these 
models  continues  to  simulate  the  transistors  without  the 
protecting  diodes  (which  are  shown  in  Figure  2.2)  in  crder 
to  demonstrate  the  inductive  voltage  transients  produced  by 
opening  one  of  the  stator  winding  circuits  with  a  transistor 
switch.  The  next  mcdel  adds  the  damper  diodes  to  the  power 
conditioner  to  fully  simulate  the  power  conditioner  final 
drive  elements.  The  final  model  takes  all  the  complications 
of  the  preceeding  ones  and  depicts  the  transient  behavior  of 
the  system  during  reversal  ("plugging")  of  the  rotating 
motor.  The  results  of  computer  runs  using  this  mcdel  are 
discussed  in  Chapter  Five.  No  attempt  has  been  made  to 
reduce  the  order  of  the  model  in  these  simulations  since  the 
desire  was  to  observe  both  the  fast  and  slow  time  constant 
effects.  The  numerical  integration  of  the  differential  equa- 
tions of  such  a  "stiff"  system  to  evaluate  the  steady  state 
behavior  can  be  somewhat  time  consuming.  Transistor  time 
constants  are  user  adjustable  as  are  saturation  and  cutoff 
parameters.  The  motcr  parameters  used  for  the  model  devel- 
opment are  typical  fcr  the  brushless  motors  now  available 
for  this  type  of  application,  but  are  not  representative  of 
any  particular  motor.  As  such,  all  motor  parameters  are 
user  adjustable  to  ccnform  to  any  desired  motor  perf crmance 
that  would  meet  design  objectives. 


19 


III.  ELECTRONIC  COMMUTATION 

A.   THREE  PHASE  BRIDGE  CIRCUITRY 

The  three-phase,  four  pole  brushless  DC  motor  is 
depicted  in  Figure  2.2.  The  transistor  and  diode  configura- 
tion is  very  similar  in  appearance  to  a  three  phase  full 
wave  bridge  rectifier  circuit  and  acts  in  a  very  similar 
fashion.  The  composite  flux  shown  in  Figure  4.6  is  a  demon- 
stration of  this  similarity.  The  composite  flux  actually 
corresponds  to  the  peaks  of  the  back  emf  voltage  waveform  in 
Figure  4.5  as  if  it  were  rectified  in  a  three  phase  full 
wave  bridge  rectifier.  This  figure  was  generated  by  taking 
the  normalized  flux  model  for  the  two  conducting  phases 
during  each  commutation  interval,  adding  them  together, 
multiplying  by  the  torgue  coefficient  and  converting  to 
webers  (the  torgue  coefficient  has  the  same  basic  units  of 
force-length/amp) .  Actually,  the  diodes  may  be  used  as  a 
full  wave  bridge  rectifier  for  some  applications  (e.g.  to 
recharge  batteries)  when  the  motor  is  driven  by  the  lead 
torgue  and  acts  as  a  generator  [Ref.  3]. 

Transistors  Q1  through  Q6  control  the  current  switching 
in  accordance  with  the  timing  diagram  in  Figure  3.1,  where  a 
logical  "1"  corresponds  to  the  transistor  being  in  the 
conduction  mode.  This  diagram  was  derived  for  counterclock- 
wise rotation  (positive  direction)  as  the  current  will  be 
positive  in  each  winding  as  the  north  pole  of  the  rotor  is 
positioned  in  its  field  (positive  flux)  and  the  current  will 
be  negative  when  the  south  pole  of  the  rotor  is  positioned 
in  its  field  (negative  flux).  A  similar  timing  diagram  can 
be  made  by  considering  the  logic  for  clockwise  rotation 
(negative  direction)  in  Table  II. 
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TRANSISTOR  SWITCHING  LOGIC 
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Figure  3.  1   Transistor  Timing  Diagram 
(Level  1  Represents  "On") 
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E.   SI1TCHING  TBANSIS10H  DYNAMICS 

Behavior  of  transistors  in  the  cutoff  and  saturation 
regions  differs  somewhat  from  the  more  typically  understood 
small  signal  models.  Some  discussion  of  switching  charac- 
teristics and  the  resulting  model  is  in  order.  This  discus- 
sion assumes  a  generic  transistor  and  is  equally  correct  for 
PNP  and  NEN  devices,  as  the  discussion  describes  primarily 
base  region  effects.  It  might  be  necessary  to  keep  in  mind 
that  electrons  are  majority  carriers  for  n-type  base  regions 
and  holes  are  majority  carriers  for  p-type  regions.  This 
discussion  of  transistor  behavior  follows  that  fcund  in 
£Eef.  4]. 

When  a  transistor  is  operating  in  the  cutoff  region,  the 
emitter  and  collector  are  reverse  biased,  base  region 
minority  carrier  concentration  is  practically  zero  at  the 
junctions,  and  each  junction  is  drawing  some  fraction  of 
leakage  current  Ico.  The  majority  carriers  which  are 
collected  at  the  junctions  are  thermally  generated  in  the 
base  region,  and  a  (small)  minority  carrier  current  will 
flow  cut  cf  the  base  terminal  to  keep  the  base  region 
neutral.  The  amount  of  minority  carrier  charge  stored  in 
the  base  region  is  negligibly  small. 

To  turn  the  transistor  on,  charge  carriers  must  be 
injected  at  the  emitter  junction  so  that  there  will  be  a 
large  charge  gradient  toward  the  collector.  This  means  that 
there  must  be  minority  carriers  stored  in  the  base  region 
and  there  must  be  a  corresponding  change  in  majority 
carriers  to  neutralize  the  stored  charge.  Thus  current 
flows  into  the  base  region  through  the  base  terminal  as  the 
storage  takes  place.  If  the  base  current  is  suddenly 
increased,  there  will  be  a  large  charge  distribution 
gradient  near  the  emitter  and  a  small  gradient  near  the 
collector,  the  difference  in  slope  being  proportional  to  the 
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ia.se  current.  Eventually  the  charge  distribution  will 
become  more  uniform  with  a  constant  gradient  of  charge 
between  the  emitter  and  collector.  This  is  the  normal 
active  transistor  region.  If  the  load  line  of  the  circuit 
is  such  that  the  transistor  is  allowed  to  pass  through  the 
active  region  and  enter  the  saturation  region,  the  collector 
current  reaches  a  maximum  value  and  cannot  increase.  The 
collector  to  base  voltage  is  then  very  small,  but  the  base 
current  can  be  increased  further  so  that  the  minority 
carrier  distribution  in  the  base  region  has  essentially  zero 
gradient  from  the  emitter  to  collector.  This  creates  a 
large  excess  of  minority  carriers  at  the  collector  junction 
which  in  turn  causes  the  collector  to  be  forward  biased. 
This  is  a  state  in  which  there  is  more  charge  at  the 
collector  than  the  collector  can  collect.  In  this  condi- 
tion, the  collector  to  emitter  voltage  is  smaller  than 
either  of  the  junction  voltages. 

The  amount  of  charge  stored  in  the  base  region  in  this 
condition  is  very  large.  Additionally,  since  the  collector 
junction  is  forward  tiased,  excess  minority  carriers  will 
cross  the  junction  from  the  base  region  and,  if  the  minority 
carrier  lifetime  is  sufficiently  long,  be  stored  in  the 
collector  region  as  well  as  in  the  base  region.  This  stored 
minority-carrier  charge  has  a  deleterious  effect  on  the 
turn-off  time  of  the  transistor.  In  order  to  turn  off  the 
transistor,  the  minority  carrier  density  at  the  collector 
junction  must  be  reduced  to  zero  for  the  collector  voltage 
to  start  to  recover.  When  turn-off  begins,  the  charge 
gradient  is  smaller  at  the  emitter  than  at  the  collector, 
resulting  in  current  flow  into  the  base  lead  and  minority 
carriers  disappearing  from  the  base  region. 

The  collector  current  remains  unchanged  until  the 
storage  period  is  over  and  collector  voltage  remains  small. 
As  the  transistor   enters  the  active  region,    the  remaining 
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charge  is  swept  out  of  the  collector  and  emitter.  The 
minority  carrier  density  at  the  collector  is  again  small, 
and  the  charge  gradient  is  decreasing  rapidly.  As  socn  as 
minority  carrier  density  reaches  zero  at  the  emitter,  it  may 
also  teccme  reverse  tiased.  This  turn-off  process  is  one  of 
the  important  limitations  on  the  switching  response  of 
transistors. 

The  turn-on  time  for  a   simple  switch  is  approximated  by 
the  following  eguaticrs:   [Ref.  4] 
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and  ^  is  the  -3dB  alpha  cutoff  frequency. 

The  effects  of  overdrive  and  reverse  drive  on  transistor 
switching  performance  are  demonstrated  in  Figure  3.2.  If 
not  overdriven,  it  takes  approximately  three  time  constants 
(Tec  or  Tee)  to  reach  saturation.  A  significant  improvement 
in  rise  time  is  realizeable  by  the  use  of  base  overdrive. 
If  a  transistor  is  overdriven,  its  effective  gain  is 
decreased,  but  this  has  the  effect  of  slightly  reducing  the 
time  constant  while,  more  significantly,  attempting  to  drive 
the  transistor  to  higher  steady  state  current — resulting  in 
a  greater  initial  slope.  The  combination  of  the  two  effects 
results  in  rise  times  on  the  order  of  .7  time  constants  to 
reach  saturation. 
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The   turn-off   time     for   a   simple    switch      is  characterized 
by  the   following   equations:      £Bef.    4] 


re=rE-e,P%7  -  7!— f'-^p^r)  (3-7) 

rc  =  Ic"exP=p-    -  ^  I   0fr(  1-  exp  ^-)  <3-8) 


and  1^  and  I~  are  currents  just  prior  to  turn-off,  and 
I8,off  is  the  base  drive  current  immediately  following  turn- 
off  initiation-  An  inspection  of  equations  3.7  and  3.8 
reveals  that  the  turn-off  time  is  affected  by  two  mecha- 
nisms, one  of  base  drive  current  and  one  of  collector  or 
emitter  current  initial  condition.  Reducing  the  initial 
condition  of  collector  or  emitter  current  or  increasing  the 
base  drive  current  results  in  a  proportional  decrease  in 
turn-cff    time. 

One  of  the  difficulties  of  overdrive  is  storage  time 
which  is  a  function  of  the  excess  base  current.  Storage 
time  for  a  simple  switch  is  described  by  equations  3.S  and 
3.  10; 

CE:     Ts  =  rso  i„  «c(v-v)  {3.io) 

zrc"  -  <%ieZq 

where       lso  =  ,  , r  (3.11) 

and  <x%  corresponds  to  or3  with  collector  and  emitter  termi- 
nals reversed,  and  a^  is  the  -3dB  alpha  cutoff  frequency 
with    collector   and   emitter   terminals   reversed.         [Ret.    1   and 

5]- 

It  is  most  desirable  to  minimize  the  storage  time  an  1 
this  can  be  accomplished  by  overdriving  during  turn-on, 
reducing   the   base   drive   to    hold   the   transistor   near 
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saturation  and  finally  drive  it  off  quickly  using  reverse 
drive.  This  can  be  accomplished  neatly  and  simply  by 
bypassing  the  base  resistor  with  a  suitably  sized  capacitor 
as  discussed  in  [Hef.  6  :p,  784]. 


C.   HCDEI  OF  SHITCHIBG  TRANSISTOR  DYNAMICS 

The  preceding  discussion  leads  into  the  modeling  of  the 
transistor  switches  as  dynamic  elements.  The  most  desirable 
transistcr  model  would  be  of  low  order  and  not  add  signifi- 
cantly tc  the  program  complexity  since  the  motor  dynamics 
are  the  prime  coccerr. 

The  initial  models  were  developed  with  the  transistor 
switches  acting  instantaneously  between  cutoff  resistance 
and  saturation  resistance.  This  was  satisfactory  as  lccg  as 
the  transients  generated  remained  small.  Using  the  variable 
step  size  integration  routines  available  with  CSMP,  it  was 
found  that  the  rapid  switching  drove  the  integration  step 
size  minutely  small  ard  computer  time  for  the  simulation  to 
reach  steady  state  became  excessively  long.  Based  upon  the 
discussion  given  in  the  preceeding  section  on  transistor 
switching,  an  exponential  rise  and  fall  was  developed,  and 
both  overdrive  and  reverse  drive  simulated  by  clippirg  the 
exponential  rise  and  decay  within  limits  of  the  saturation 
resistance  and  cutoff  resistance.  This  not  only  provided  a 
good  simulation  of  the  transistor  switching  with  user  vari- 
able parameters,  but  also  eased  the  computer  time  problem 
somewhat.  The  resulting  switching  characteristics  are  shown 
in  Figure  3.3,  (with  exaggerated  time  constants  for  clarity) 
and  compare  favorably  with  those  shown  in  Figure  3.2  (a)  and 
(b) .  Figure  2.2  suggests  that  the  transistors  are  bipolar 
NPN  devices.    This  was  done   for  pictorial   reference  only. 
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Figure   3.2 


Switching  Characteristics,  (aj  Overdrive  only, 
(b)  Overdrive  and  Reverse  Drive. 
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TABLE 

III 

Some  Power 

Transistor  Specifications 

Device 

Ic  (flax) 
Amps 

H  (sat) 

Ohms 

ts 
s 

Diss 
Watts 

Vbr  (sus) 
Volts 

MTM45N15 
(BCS) 

45* 

.065? 

250 

150 

MTM40N18 
(MOS) 

40* 

.08% 

250 

180 

MJ16016 
(Eipclar) 

20 

2.2 

250 

450 

HJ 15011 
(Eipclar) 

10 

200 

250 

2U3430 
(Eipclar) 

5 

.2 

4.0 

150 

100 

2N2777 
(Eipclar) 

30 

.06 

14** 

200 

150 

2N6250 
(Eipclar) 

15 

3.5 

175 

275 

2K6339 
(Eipclar) 

25 

1.0 

200 

120 

2N1820 
(Eipclar) 

15 

.1 

25** 

250 

250 

STC2231 
(Eipclar) 

30 

.08 

6# 

200 

200 

♦ID 
3E  (on) 

**ts+tf 

#td+tr 

J 
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These  switching  transistors  may  be  a  pair  of  complementary 
symmetry  Class  B  amplifier  output  devices  to  take  advantage 
of  emitter-follower  characteristics  or  might  be  MOS  devices 
to  realize  high  trarsistor  power  efficiency.  The  model  as 
presented  is  accurate  for  bipolar  transistor  dynamics  tut 
because  switching  time  constant  and  on  and  off  resistance 
are  the  cnly  parameters  for  the  model,  it  should  be  exten- 
dable to  include  MOS  devices  as  well.  The  parameters  used 
for  the  simulation  prcgram  were  a  switching  time  constant  of 
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1.5  microseconds  and  a  saturation  resistance  of  .05  ohms. 
These  values  were. chosen  as  more  or  less  ideal  values,  the 
combination  of  which  would  probably  not  be  realized  in  any 
one  device.  In  Table  III  some  transistor  specifications  are 
listed  for  devices  which  might  be  desirable  to  use  for  this 
motor  controller  application.  The  fast  switching  times  for 
the  bipolar  transistors  are  their  salient  feature,  while  the 
low  "on"  resistance  of  the  MOS  devices  makes  them  desirable. 


TRANSISTOR  SATURATION  MODEL 
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Figure  3.3   Transistor  Model  Saturation  Characteristics 

The  model  developed  to  this  stage  has  ignored  the  effect 
of  peak  reverse  voltage  across  the  transistors  that  occurs 
at  the  moment  of  ccnmutation  due  to  the  opening  of  the 
inductive   circuit.    To   provide   the  required   protection, 


29 


damper  diodes  are  used  to  clamp  the  spike  voltage  tc  the 
positive  cr  negative  supply.  The  diodes  not  only  provide 
overvcltage  protection  for  the  transistors,  but  must  handle 
nearly  the  full  value  of  steady  state  current  for  a  brief 
period  of  time  as  well.  The  effect  of  the  diode  clamping 
can  fc€  seen  immediately  by  comparing  the  figures  depicting 
motor  current  with  and  without  the  diodes  in  the  simulation. 
(See  Figures  5.1  and  5. 2) Note  that  the  current  dip  is  mark- 
edly reduced  when  the  diode  is  in  the  circuit.  The  diodes 
provide  a  commutation  of  the  transistor  current  and  mairtain 
current  flow  in  the  motor  during  the  switching  of  the  drive 
current  between  two  phases.  additional  information 
concerning  diode  commutation  was  found  in  [Eef.  9  :p.  371] 
concerning  the  so-called  "Jones"  commutation  circuit  for  a 
EC  mctci  controller. 

The  diode  model  is  very  simple  and  features  instanta- 
neous switching  between  saturation  and  cutoff  resistance 
identical  in  value  to  that  of  the  transistors.  This 
approach  did  not  add  any  appreciable  complexity  to  the  inte- 
gration problem  because  the  diodes  pick  up  the  current  as 
the  transistor  turns  off,  smoothing  out  the  transient 
current  instead  of  generating  another  current  transient. 
The  switch  of  the  current  to  the  diode  is  triggered  by  the 
voltage  spike  caused  by  the  transistor  turning  off,  similar 
to  what  would  be  expected  in  the  actual  circuit.  The  lack  of 
a  diode  voltage  drop  before  the  knee  of  the  diode  conduction 
curve  is  the  only  addition  that  might  be  made  to  the  model, 
but  was  deemed  unnecessary  as  the  voltage  drop  across  the 
diode  at  full  meter  current  would  still  be  less  than  one 
volt  and  of  the  same  crder  of  magnitude. 
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IV.  MAGNETIC  FLOX  MODEL 

A.   STMMETBIC  FIFTH  BABHONIC  FIOX 

Thomas  [Bef.  1]  derived  a  combination  fundamental  and 
fifth  harmonic  reverse  voltage  waveform  from  empirical  meas- 
urements. This  model  has  been  used  in  this  thesis  with 
modifications.  The  phase  relationship  of  the  fifth  harmonic 
to  the  fundamental  has  been  altered  to  provide  the  same 
shape  voltage  waveform  across  two  windings  instead  of  across 
only  one.  This  agrees  with  the  empirical  data  and  also  gives 
rise  to  a  new  magnetic  flux  model,  the  shape  of  which  agrees 
in  form  with  traditicral  DC  machines  [Bef.  8].  The  general 
shape  of  the  flux-density  distribution  of  a  DC  motor  is 
shown  in  Figure  4.1. 
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Figure  4.1    Flux  Density  Distribution  due  to  Botor  Field 
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The  dip  in  the  center  of  the  flux-density  distr ituticr  can 
test  re  explained  by  considering  Figures  4.2  and  4. 3  which 
are  flux  maps  of  the  magnetic  field  due  to  the  rotcr  magnets 
and  stater  current  respectively  and  confirm  the  distribution 
shown  in  Figure  4.1  for  either  stator  or  rotor  considered. 
In  Figure  4.2  and  Figure  4.3,  separation  of  flux  lines  indi- 
cates degree  of  field  strength,  with  greater  separation 
indicating  less  field  strength. 


Figure  4-2    Flux  Bap  with  only  Rotor  Field  Considered 
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Figure  4.3   Flux  Hap  with  only  Stator  Field  Considered 


The  model  flax-density  distribution  is  depicted  in 
Figure  4.4  showing  the  phase  relationship  between  the  three 
windings.  In  Figure  4.5  the  phase  fA'  flux  is  shown  alcng 
with  th€  resulting  generated  voltage  between  terminals  'A' 
and  *B'.  The  dotted  lines  in  each  figure  represent  the 
commutation  periods  during  which  each  winding  is  connected 
to  the  positive  supply,  the  negative  supply,  or  effectively 
disconnected  from  supply  by  the  switching  power  transistors. 
The  commutation  effects  shown  by  the  dotted  lines  in  these 
figures  represent  motcr  currents  for  counterclockwise  (posi- 
tive) rotation  as  the  product  of  the  flux  and  the  motor 
current  is  seen  to  always  be  positive.  Since  the  motor 
torque  is  derived  from  the  product  of  the  magnetic  flux 
density  distributions  and  the  current  flowing  in  each  phase 
winding,  this  results  in  a  positive  motor  torque  and  rota- 
tion in  the  positive  direction.  The  composite  flux  density 
distrituticn  for  the  rotating  motor  is  displayed  in  Figure 
4.6,  and  is  used  to  derive  the  normalizing  constant  for  the 
flux  density  by  maintaining  the  average  value  for  the  torque 
coefficient  at  the  constant  value  for  the  simple  models. 
This  composite  flux  density  distribution  is  seen  to  have  an 
average  value  equal  to  the  torque  constant  of  the  simple 
models,  when  the  conversion  from  webers  to  inch-ounces  per 
amp  is  applied: 

(in-cz/amp)  x  (.0C707)  =  weters  (4.1) 

This  figure  also  provides  a  quick  check  on  the  back  eraf 
constant  since  the  numerical  value  of  both  constants  must  be 
identical  in  HKS  units. 

To  demonstrate  clockwise  rotation,  the  flux  model 
remains  fixed  with  respect  to  the  shaft  position,  but  the 
sign  of  the  currents  must  be  reversed.  This  amounts  to 
reversing   the  polarity   of   the   commutation  scheme   for   a 
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Figure    4.4        Model   Flux— Phase  Relationship 
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MAGNETIC  FLUX  MODEL 
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Figure  4-5    Model  Flux  and  Back  EMF  Voltage 
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MAGNETIC  FLUX  MODEL 

COMPOSITE  FLUX  VARIATION 
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Figure  4.6    Hodel  Composite  Flux 

counterclockwise  rotating  motor.  As  can  be  seen  frca  Table 
II  this  is  the  same  sequence  of  commutation  as  for  counter- 
clockwise rotation  advanced  by  ninety  degrees,  which  should 
be  expected  from  the  rotor  geometry  since  the  south  pcle  of 
the  rotor  is  displaced  ninety  degrees  from  the  north  pcle  of 
the  rctor. 


E.   FIFTH  HARMONIC  FI0X  WITH  DEIVE  CURRENT  INTERACTION 

The  interaction  cf  the  two  magnetic  fields,  one  due  to 
the  rctcr  magnets  and  the  other  due  to  the  stator  current 
causes   distortion   of   the   air   gap   flux   and   hence   the 
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generated  voltage  waveform.  Reference  eight  indicates  that 
such  an  interaction  is  expected,  but  as  to  what  degree  the 
flux  density  distribution  and  the  reverse  voltage  waveform 
would  be  distorted  is  unknown.  The  expected  effect  of  the 
two  magnetic  fields  interacting  is  to  skew  the  symmetrical 
distribution  as  is  shewn  in  Figure  4.7  It  is  realized  that 
the  profile  of  the  flux  model  of  Figure  4.5  and  the  total 
flux  density  distribution  as  shown  in  Figure  4.7  do  rot 
agree,  however  the  flux  model  shows  a  high  degree  of  corre- 
lation to  empirically  measured  experimental  test  data  for 
this  type  cf  motor  [Bef.  11]. 
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Figure  4.7   Flux  Density  Distribution 
Demonstrating  Interaction. 


The  skewed  shape  cf  the  resultant  flux-density  distribu- 
tion due  to  stator  current  interaction  can  be  simulated  by 
slipping  the  phase   relationship  of  the  fifth   harmonic  with 
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respect  to  the  fundamental.  If  the  phase  shift  is  simply  a 
constant,  then  two  windings  will  produce  essentially  the 
same  shape  of  hack  EHF  voltage  as  the  previous  model,  tut 
combination  with  the  other  winding  will  not  be  similar  or 
symmetric.  Since  current,  as  modeled,  essentially  flows  in 
two  windings  at  a  time,  the  distortion  must  be  a  function  of 
the  current  and  the  third  winding  flux  must  remain  undis- 
torted  as  long  as  that  winding  has  no  current  flowing. 
This  thinking  led  to  the  model  for  interactive  flux  depicted 
in  Figure  4.8.  The  phase  relationship  between  fundamental 
and  harmonic  is  identical  to  the  non-interactive  model  if 
the  drive  current  is  zero.  In  the  presence  of  drive 
current,  the  fifth  harmonic  is  shifted  smoothly  and  continu- 
ously until  at  maximum  current,  the  phase  shift  is  thirty 
electrical  degrees.  Even  with  this  visible  warping  cf  the 
flux,  the  back  emf  voltage  waveform  is  recognizeatle  as  a 
fundamental  and  a  fifth  harmonic,  and  this  effect  might  be 
smaller  than  depicted,  and  of  little  consequence  tc  the 
motor  designer.  Additionally,  no  model  of  magnetic  satura- 
tion has  been  implemented,  which  would  have  the  effect  of 
flattening  the  flux  peaks  and  giving  a  more  constant  maximum 
flux  and  hence  torgue,  resulting  in  smoother  motoring.  Such 
effects  are  expected  in  an  actual  motor,  but  at  what  value 
of  current  the  saturation  should  begin  will  have  tc  be 
determined  by  measurement  of  a  particular  motor. 
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Figure  4.8    Model  Flux  and  Back  EMF  Voltage 
Demonstrating  Drive  Current  Interaction 
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1.  HOD EL  PEBFOBMANCE 

The  performance  cf  the  first  three  models  is  net  inter- 
esting to  the  motor  designer  as  the  simplifications  involved 
in  these  simulations  reduce  the  resulting  behavior  to  text- 
book DC  motor  dynamics.  As  a  development  tool  for  other 
simulations,  the  structure  may  be  of  some  interest.  These 
simulations  are  included  in  Appendix  C  for  this  reason.  The 
transistor  switching  of  the  fourth  and  all  subsequent  models 
is  of  interest  in  that  the  transistor  impedance  affects  the 
(L/R)  time  constants  of  the  motor  current  and  is  an  impor- 
tant feature  in  understanding  the  behavior  of  the  simula- 
tions. Because  of  the  varying  time  constants,  decay  of 
current  in  one  phase  and  buildup  of  current  in  another  phase 
at  commutation  is  net  compensatory.  The  current  decays 
faster  in  the  turn-off  winding  than  the  corresponding 
buildup  in  the  turn-cn  winding,  while  the  current  in  the 
third  winding  is  ostensibly  constant.  The  third  winding 
current  cannot  be  constant  because  all  three  currents  oust 
sum  to  identically  zero.  The  result  is  the  sharp  dip  in  the 
center  of  the  current  waveform  of  the  unswitched  winding. 
Excellent  correlation  of  this  effect  is  seen  in  the  experi- 
mental data  as  shown  in  Figure  (12)  of  [Ref.  10]  and 
[Bef.  11]-  In  simulations  which  include  a  diode  model,  the 
dynamic  equations  that  determine  the  (L/E)  time  constant  are 
expanded  to  include  the  diode  model.  The  inclusion  cf  the 
diode  impedance  in  the  dynamic  time  constant  results  in  a 
"bracing  up"  of  the  current  waveform  of  the  unswitched 
winding  yeilding  a  higher  average  motor  current.  This  can  be 
observed  by  comparing  Figures  5.1  and  5.2.  Models  with 
simulated  diodes  demonstrate  somewhat  more  torque  and  speed 
with  slightly  less  speed  ripple  than  those  without  diodes. 
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The  final  versions,  in  which  the  transistor  and  diode 
models  and  sinusoidal  back  emf  voltage  are  incorporated  are 
used  to  characterize  an  actual  motor  and  power  conditioner. 
The  performance  curves  provided  by  Thomas  [Ref.  1]  were  used 
for  comparison.  As  a  result,  all  simulations  presented  here 
used  a  net  supply  voltage  of  thirty  volts,  or  plus  and  minus 
fifteen  volts  applied  to  the  transistors.  As  stated  earlier 
the  object  of  this  tiiesis  was  to  provide  a  tool  to  evaluate 
motor  and  switching  transistor  performance  for  various 
design  parameter  values  and  operating  conditions.  Seme 
deviation  from  these  performance  curves  was  anticipated,  tut 
overall  characteristics  of  the  torque,  current,  and  power  vs 
speed  curves  remain  of  the  same  form. 

At  the  time  of  this  writing,  provisions  are  being  made 
to  gather  mere  comprehensive  test  data  from  a  typical  motor. 
It  is  well  understood  that  only  a  limited  number  of  the 
model  waveforms  may  be  confirmed  by  measurement.  The  brush- 
less  DC  notor  is  a  three  wire  device,  and  the  voltage  across 
any  twe  windings  may  be  conveniently  measured.  The  simula- 
tion, however,  calculates  and  graphs  not  only  the  voltage 
across  each  pair  of  windings  but  in  addition,  the  voltage 
across  each  single  winding.  The  individual  currents  IA,  IB, 
and  IC  may  also  be  measured  with,  for  example,  an  inductive 
coupler  or  by  measuring  the  voltage  drop  across  a  resistor. 
These  measured  curves  can  be  readily  compared  with  model 
results  and  model  parameters  adjusted  to  bring  about  good 
agreement. 
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Mctor  reversal  at  speed  can  be  simulated  under  varying 
load  conditions.  In  the  simulation  the  load  torque  opposes 
motor  rotation,  and  so  helps  the  motor  to  slow  to  a  stop, 
then  cpposes  the  reversal.   If   a  viscous  load  was  simulated 

(as  in  driving  a  dynamometer) ,  slowing  to  a  stop  would  take 
longer  because  the  lead  would  decrease  with  speed,  and  the 
rise  to  steady  state  in  the  reverse  direction  wculd  be 
quicker  due  to  the  reduced  load  at  low  speed.  In  general, 
this  type  of  reversal  is  characterized  by  large  trarsient 
motor  currents  since  the  back  emf  voltage  is  initially  of 
polarity  that  aids  current  flow  instead  of  inhibiting  it. 
Ihe  reverse  torque  generated  by  this  current  will  quickly 
slow  the  motor  to  a  stop  where  there  is  no  back  emf  gener- 
ated, and  as  the  rotation  reverses  the  generated  voltage 
changes  polarity  and  once  again  opposes  current  flew.  The 
reversal  effects  can  be  seen  in  the  current  waveforms  of 
Figures  5.3,   5.4,   5.5  and   5.6  which  depict  motor  reversal 

(from  counterclockwise  rotation  to  clockwise  rotation)  at 
time  =  . C2  seconds. 
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In  these  figures,  a  demand  for  reversal  is  initiated  at 
.02  seconds.  Since  the  rotational  inertia  of  the  rotor  keeps 
it  rotating  in  the  positive  direction  for  a  short  time  after 
reverse  current  begins  to  flow,  a  commutation  point  is 
reached  and  the  reversal  current  is  switched  between  phases 
A  and  C.  As  the  rotor  stops  and  the  motor  reverses  under  the 
influence  of  negative  torque,  the  commutation  point  is  again 
reached  hut  from  the  opposite  direction  and  the  current  is 
switched  from  phase  C  to  A.  There  is  a  large  current  tran- 
sient in  the  figure  at  this  point  which  is  also  seen  in 
Figure  5.7  as  a  transient  positive  torque,  probably  due  to  a 
combination  of  the  irstantanecus  nature  of  the  diode  model 
and  the  transistor  time  constants.  The  reality  of  such  a 
current  and  torque  spike  is  questionable,  but  if  the  dynamic 
nature  of  the  transistors  and  diodes  used  in  the  power 
conditioner  are  not  similar,  such  undesirable  transients  may 
be  generated,  even  to  the  point  of  punch  through  of  a  junc- 
tion. Ihe  effect  of  the  torque  spike  is  seen  in  Figure  5.6 
as  a  slight  interruption  of  the  otherwise  smooth  transition 
between  forward  rotation  and  reverse  rotation. 

The  simulation  can  be  a  powerful  analytic  and  design 
tool.  For  example,  the  addition  of  damper  diodes  was  seen 
to  maintain  current  flow  during  transistor  switching  (demon- 
strated in  Figures  5.1  and  5.2)  which  has  the  effect  of 
maintaining  a  higher  average  torque  with  diodes  than  without 
as  has  been  previously  discussed.  It  follows  that  the 
faster  the  diodes  operate,  the  better  the  motor  performance 
which  can  be  realized  in  this  context. 

The  availability  of  the  center  of  the  three  windings  for 
use  as  a  ground  terminal  may  prove  to  be  advantageous  in 
some  applications.  The  current  return  path  to  the  power 
supply  provided  by  the  center  connection  would  not  be  opened 
by  the  switching  transistors,  and  the  current  flow  *ould  not 
be  interrupted  in  the  unswitched  winding  during  commutaticn, 
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resulting  in  smoother  rotation.  The  design  tradeoff  to  this 
would  be  that  only  half  of  the  supply  voltage  would  be 
available  across  any  one  winding  with  the  center  grounded. 
If  the  center  terminal  is  left  to  float,  the  voltage  across 
two  windings  is  the  full  supply  voltage,  but  the  voltage 
applied  across  a  single  winding  and  the  current  through  that 
winding  is  a  function  of  the  back  emf  voltage  generated  in 
both  windings.  Since  the  back  emf  voltage  is  varying  siru- 
soidally  in  each  winding  and  one  winding  is  120  electrical 
degrees  cut  of  phase  with  respect  to  each  of  the  other  two 
windings,  a  larger  proportion  of  the  supply  voltage  may  be 
applied  to  one  of  the  conducting  coils  than  the  other.  The 
current  through  the  cne  coil,  being  equal  to  the  sum  of  the 
currents  in  the  other  two  coils,  is  at  a  higher  value  when 
the  air  gap  flux  is  greatest  than  if  the  center  connection 
were  grounded,  and  will  provide  more  motor  torque.  The 
trajectory  of  the  center  connection  is  shown  in  Figure  5.8 
and  can  be  compared  tc  the  single  phase  flux  characteristic 
in  Figure  4.4  and  the  composite  flux  variation  in  Figure  4.6 
to  see  that  the  voltage  aids  the  current  flow  in  each 
winding  as  it  comes  under  the  peak  of  the  magnetic  flux,  as 
the  voltage  at  the  center  connection  is  negative  at  during 
the  positive  flux  peaks,  and  positive  during  the  negative 
flux  peaks,  aiding  current  flow  when  compared  to  a  grounded 
center  terminal.  The  relative  scale  on  the  ordinate  of  this 
figure  indicates  that  the  effect  is  proportional  tc  motor 
speed. 

The  fluctuation  of  motor  speed  demonstrated  by  this 
model  has  been  of  the  order  of  ten  percent  of  peak  motor 
speed,  somewhat  more  than  in  Thomas'  simulations  [Bef.  1]. 
This  can  be  attributed  to  the  combination  of  the  transistor 
switching  time  constants,  the  current  rise  time  constant 
(L/R)  and  extremely  small  rotor  inertia.  If  the  load 
inertia  was  of  the  order  of   magnitude  of  the  rotor  inertia, 
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CENTER  POINT  TRAJECTORY 
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Figure  5.8    Trajectory  of  Center  Connection 

the  variation  in  speed  would  te  reduced  to  the  order  of  five 
percent  due  to  tie  increased  time  constant,  («VE)  . 
Additionally,  the  speed  ripple  decreases  significantly  if 
the  applied  voltage  is  increased  with  no  increase  in  load, 
flith  a  load  torque  of  64  inch  ounces  and  a  power  supjiy 
voltage  of  sixty  vclts,  speed  ripple  is  hardly  noticeable  at 
the  steady  state  speed  of  approximately  4800  rpm  due  to  the 
increased  angular  momentum. 
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with  minor  changes  to  the  program,  different  magnetic 
flux  shapes  can  he  simulated  and  motor  performance  can  be 
compared.  For  example,  if  the  stator  was  designed  with 
interpoles,  the  flux  would  take  on  a  flatter  shape  across 
the  top  of  the  flux  curve,  which  would  in  turn  give 
different  characteristics  to  the  back  emf  voltage,  which 
could  be  empirically  measured.  If  the  torque  coefficient  of 
the  motor  was  known,  (it  can  be  measured)  the  flux  distribu- 
tion could  be  deduced  as  it  has  been  for  these  simulations 
and  the  complete  motcr  performance  curves  predicted.  The 
simulations  give  a  good  representation  of  transient  and 
steady  state  motor  and  transistor  currents  [Ref.  10],  from 
which  voltage  and  current  stresses  applied  to  the  transistor 
junctions  can  be  determined.  The  results  of  a  simulation 
for  a  desired  steady  state  load  may  be  analyzed  to  size  the 
transistors  and  diodes  required  for  a  specific  application. 
Given  the  heat  transfer  characteristics  of  a  motor  and/or 
heat  sink  assembly,  the  thermal  rise  of  the  power  condi- 
tioner and  the  motor  stator  may  be  predicted  by  implementing 
a  simple  heat  transfer  model. 

from  the  literature  [Bef.  3]  it  was  noted  that  when  some 
thirty  degrees  of  ccirmutation  advance  was  used  to  obtain 
maximum  performance  from  the  motors  being  tested,  an 
increase  of  2303  over  rated  horsepower  was  realized  for 
short  periods.  This  simulation  program  has  the  capability 
of  shifting  the  ccnmutation  angles  to  realize  extended 
performance.  This  is  left  for  future  development  of  motor 
performance  studies. 
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71.  CONCLUSION  AN  I)  RECOMMENDATIONS 

A.  SEMMARI  AND  CONCLUSION 

The  development  cf  a  detailed  transient  model  for  simu- 
lation cf  the  instantaneous  performance  of  a  small  fcur  pcle 
hrushless  DC  motor  supplied  by  a  transistorized  power  condi- 
tioner has  been  presented.  This  model  was  part  cf  a  compre- 
hensive development  cf  a  fin  control  actuating  system  model 
for  the  simulation  cf  the  dynamic  performance  cf  a  cruise 
missile  fin  controller. 

The  results  of  the  numerical  simulation  of  the  instanta- 
neous meter  and  power  conditioner  voltages  and  currents 
demonstrated  a  great  deal  of  correlation  to  reported  experi- 
mental test  data  durirg  steady  state  operation.  The  dynamic 
reversal  of  the  model  was  performed.  The  corresponding 
dynamic  currents  predicted  by  the  model  were  analyzed  and 
large  transients  were  demonstrated  as  would  be  expected  in 
an  actual  motor.  The  conclusion  that  the  model  is  a  good 
representation  of  an  actual  motor  and  power  conditioner 
drive  transistors  is  supported  by  all  currently  available 
motor  data  and  analytical  extensions  thereof. 

B.  RECOMMENDATIONS 

This  model  will  be  further  improved  by  the  addition  of 
delay  between  the  Hall  sensors  and  the  commutation  algorithm 
to  more  realistically  simulate  the  electronic  logic  cf  the 
power  conditioner.  A  routine  for  the  calculation  cf  tran- 
sistor and  diode  power  dissipation  and  possibly  motor  stator 
heating  will  be  included.  This  will  allow  maximum  perform- 
ance to  be  extracted  from  the  motor  and  power  conditioner 
for  short  periods  of  time  without  excessive  heating. 
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Additional  work  en  the  dynamic  magnetic  flux  profile, 
particularly  in  the  area  of  magnetic  saturation  of  the 
stator  iron  needs  tc  be  investigated  and  implemented  to 
accurately  simulate  high  power  operation. 
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APPENDIX  A 
SPACI  DEPENDENT  FLOX  MODEL 


Because  of  the  length  and   complexity  of  the  CSMP  model, 

a  simpler   model   was   designed   using   Digital   Simulation 

language   (DSL) .    The   advantages  to  DSL  include   complete 

compatibility   with  CSMF   and   a   more  or   less   interactive 

running  capability  which  the  NFS  installation  of  CSMF  lacks. 

If  the  complete  model  had  been  written  in  DSL,  the  advantage 

of  quick  turnaround   of  the  simulation  would   have  been  lest 

due  to  the  running  time  of   the  program.   The  NPS  version  of 

DSL  dees  not  have  tie   STIFF  integration  routine   used  with 

CSMP   and   the   Bunge-Kutta    integration   routine   in   both 

languages  had  some   computation  problems  with  the   fast  time 

constants  of  the  power  transistors.    The  routines  contained 

in  this  appendix  were  developed  using  [Ref.  14  and  15]. 

TITLE  THESIS  PBOJECT 

TITLE  BECSH1ESS  DC  MCTCR  POWEE  CONDITIONER 

STCEAG  Z  1(500)  ,Z2(5  0  0),Z3  (500)  ,Z4  (5  0  0)  ,  Z5  (50  0)  ,  Z6  (5  00) 

INITIAL 

CONST  PI  =  3.14159#  VIF=15.,  VIB  =  -15. 

********************************* ****** 

DYNAMIC 
*************************************** 

THETA  =  RAMP  (0.0) 

THETAD  =  THETA*180./PI 

EEflFA  =  {3.*SIN  (2.*THETA  +  (11*PI/6)  ) 


U.59*SIN  M0.*THETA+(1*PI/6)  )) 

EEMFE  =  (3.  *SIN  (2.*THETA+ (7*PI/6 
1+.59*SIN (10.*TEETA+(5*PI/6) [f 


)) 


3.  *S:   (2.*THETA+ (3*PI/6) 
-59  +  SIS J10.*THETA+  (9*PI/6) ) ) 

ELIFF=  BEMFA  -  EEMFB 


45 
46 


BEEFC  =(3 
1  + 

A 

3CDIFF=  BZMFB  -  EEMFC 

CADIFF=  BEMFC  -  BEMFA 

IF  (THETAD    .GE.  180.)    GO    TC    45 

THCON    =    THETAD 

GO    TO    46 

THCCN  =  THETAD-  180. 

CCKTINUE 

0..AND.  THCON 


50 


IF  (THCON 

IF  (THCON 
IF  (THCON 
IF  (THCON 
IF  (THCON 
IF  (THCON 

sal  =  o. 

Sw2  =  0. 
SW3  =  0. 


GE, 
GE, 
GE, 
GE, 
GE.120.  .AND. 


.LT. 


30. .AND.  TECON  .  LT. 

6C..AND.  THCON  . LT. 

90. .AND.  THCON  .LT.120, 

ECON  .LT. 150. 


30 
60 
90 


GE. 150. .AND.  THCON  -LT.180 


GO  TO  50 

GO  TO  51 

GO  TO  52 

GO  TC  53 

GO  TO  5  4 

GO  TO  55 


56 


sum 

1- 

SW5    = 

1. 

SW6    = 

0. 

BEMET 

=    BEMEC    -    EEMEB 

7N1    = 

7IE    -    7EMFC 

7N2    = 

7IB    +    7EMPB 

GO    TO 

60 

51 

SW1    = 

1. 

Sfl2    = 

0. 

SH3    = 

0. 

swu  = 

1. 

SW5    = 

0. 

SW€    = 

0. 

BEMFI 

=    BEMFA    -    EEMEB 

7N1    = 

7IF    -    7EMEA 

7N2    = 

7IB    +    7EMEB 

GO    TO 

60 

52 

SH1    = 

1. 

SW2    = 

0. 

SW3    = 

0. 

SW4    = 

0. 

SW5    = 

0- 

ST76    = 

1. 

BEMET 

=    BEMFA    -    BEMFC 

7N1    = 

7IF    -    7EMEA 

7N2    = 

7IB    +    7EMIC 

GO    TO 

60 

53 

SS1    = 

0. 

Sfl2    = 

0. 

SX3    = 

1. 

ST"*    = 

0. 

SW5    = 

0. 

SW6    = 

1. 

BEHFT 

=    EEMFE    -    EEMFC 

7N1    = 

7IF    -    7EMIB 

7N2    = 

7IB    +    7EBFC 

GO    TO 

60 

54 

SW1    = 

0. 

SW2   = 

1. 

SU3   = 

1. 

SE4    = 

0, 

SW5    = 

0. 

SW6    = 

0. 

BEEFT 

=    BEMFB    -    BEMFA 

7N1    = 

7IF    -    7EMPB 

7N2    = 

7IB    +    7EMEA 

GO    TO 

60 
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SW1    = 

0. 

SW2    = 

1. 

Sfl3    = 

0. 

SW4    = 

0. 

SW5    = 

1. 

Stf€    = 

0. 

BEMFT 

=    EEMFC    -    BEMFA 

7N1    = 

7IF    -    7EMEC 

7N2    = 

7IB    +    YEHIA 

60 

CONTINUE 

IA    =    ' 

4.* 

[SW1-SW2) 

IB   = 

^.*  | 

Sfl3-Sfl4) 

IC   = 

4.* 

SW5-SK6) 

************* 

e* ******* **************** 

SAMP  IE 

*************************************** 

*      Gil    TEE 

NUMBERS 

Z1  I 

I): 

=    THETAD 

Z2  I 

I) 

=    BEMFA 

Z3  | 

I): 

=    BEMFB 

Z4  I 

Ij 

=    BEMFC 

Z5 

I); 

=    I? 

I 

57 


20 

CONTEI 

END 

STOP 


Z6  (I)  = 

Z7  1  = 
Z8jl[= 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
FOEMAT 
1=1+1 
PINTIM=6.28, 


L£LT=.0  1 


DELS=.01 
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>gr 


ra 

Z 

z 


m: 
2(629 


Begin  DISSPLA  pro 
HEAL  Z1 (529 
BEAI  z7(629 
1=0 

1  1=1+1 

BEAD  (30,5,END  =  3)  Z1 
BEAD(31,5,END  = 
READ  ('32, 5f  END= 
BEAD i 33,5, END= 
READ  '34,5, END  = 
BEAD(35,5,END= 
READ  i'36,5,END= 
BEAD(37,5.END  = 

5      FCEMAT(F20.6) 


(629)  ,Z3  (629)  ,  Z4  (629)  ,Z5  (62  9)  .  Z6  (6  29) 
8(629)  ,Z9(629)  ,Z10(629)  ,Z11  (6  29) 


CCNVEESICN 


1 


FBOM  IN 
:  Z2  (I)  * 

t' 


Z3 
Z4 
Z5 
Z6 
Z7 


I 


-CZ/AMP  TO 

2.54*9.8/ 

2.54*9.8/ 

2.54*9.8/ 

6C. 

60. 

60. 


WEBEBS: 
100.*16.*2.2 
100.*16.*2.2 
100.*16.*2.2 


♦16.82/5. 1393 
♦16.82/5.1393 
♦16.82/5. 1393 


CONTINUE 
CAII  TEK618 

SHEEPA  (' 
PAGE  (8.5 
NOBRDE 
PHYSOR(1 


CALL 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 
CAII 

CAII 
S1C£ 
ENE 


( 


(' 


AREA2D 

COMPLX 

BEADIN 

HEADIN 

XNAME 

YNAME 

XINTAX 

GEAF 

CURVE 

DOT 

CURVE 

RESET 

ENDGR 

PHYSOR('1. 

AREA2D  ( 

COMPIX 

XNAME  ('P 

YNAME  (»F 

XINTAX 

GEAF  (0.0 

CURVE  (Z1 

DOT 

CUEVE 

BESET 

ENDGE 

PHYSOE('1 


(0.0 
(Z1 

(Z1 

(•DO 

JQ1 


(Z1 

(f  DO 

m 


AEEA2D 

XINTAX 

XNAME 

YNAME 

GEAF  ( 

CUEVE 

DOT 

CUEVE 

EESET  ('"DO 

ENDPI  (0) 

DONEPI 


(5 

'E 
i 

.0 
(Z1 


MACMIILA',  •A1  ,3) 
,11.0) 

,  5.5) 

5.0,2.0) 

MAGNETIC    FIUX    MODELS ', 100, 1. 5,2) 
EEASE   A    FLUX$', 100, 1., 2) 
EASE    B    FLUX$» ,100) 
$«,100) 

,30., 36  0., -0.10, ■ SCALE', 0. 10) 
,Z2, 629,0) 

,Z5, 629,0) 

I1) 

.    3.25) 
5.0,2.0) 

HASE    C    FLUX$',100) 
IUX,    WEBEESS' , 100) 

,30.. 36  0. .-0.10, 'SCALE ',0.10) 
,Z3,o29,0f 

,Z6, 629,0) 
T») 

10*2.0) 

CTATION    ANGIE,     DE3EEES3* , 100) 

$'.100) 
,30.-36  0. ,-0. 10, 'SCALE' ,0.10) 
,Z4,S29    ,6) 

,Z7, 629,0) 
I') 
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